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Abstract 

It is of considerable importance to have a numerical method for solving supersymmetric 
theories that can support a non-zero central charge. The central charge in supersym- 
metric theories is in general a boundary integral and therefore vanishes when one uses 
periodic boundary conditions. One is therefore prevented from studying BPS states 
in the standard supersymmetric formulation of DLCQ (SDLCQ). We present a novel 
formulation of SDLCQ where the fields satisfy anti-periodic boundary conditions. The 
Hamiltonian is written as the anti-commutator of two charges, as in SDLCQ. The anti- 
periodic SDLCQ we consider breaks supersymmetry at finite resolution, but requires 
no renormalization and becomes supersymmetric in the continuum limit. In principle, 
this method could be used to study BPS states. However, we find its convergence to 
be disappointingly slow. 



1 Introduction 



Supersymmetric gauge theories in low dimensions have been shown to be related to non- 
perturbative objects in M/string theory Q, and are therefore of particular interest nowadays. 
More dramatically, there is a growing body of evidence suggesting that gauged matrix models 
in + 1 and 1 + 1 dimensions may offer a non-perturbative formulation of string theory 0, [| . 
There is also a suggestion that large iV gauge theories in various dimensions may be related 
to theories with gravity [|, 0, [J. 

It is therefore interesting to study directly the non-perturbative properties of a class of 
supersymmetric matrix models at finite and large N c , where N c is the number of gauge 
colors. Supersymmetric Discrete Light Cone Quantization (SDLCQ) [0 is a unique non- 
perturbative numerical approximation that is manifestly supersymmetric at each stage of 
the calculation. In simplest terms SDLCQ is the approximation that arises when a theory 
is confined to a box of length L in the spatial light-cone directions. This leads to a discrete 
Fock space basis and the supercharges are finite-dimensional matrices in this basis. This 
combination of the DLCQ method and supersymmetry is well defined and has proven to 
be a powerful tool which has allowed to solve a large class of problem that have not been 
solved previously [0, [D], [L4], pp . It appears that supersymmetric theories are completely 
well defined, when formulated on a compact space |13|, |15|| . All of the models that have been 
addressed so far have had sufficient supersymmetry to make them completely finite so that 
no renormalization was necessary. It is important to note that numerical approximations 
that do break supersymmetry will still be faced with a renormalization problem. 

In the light-cone formulation the supercharges and Q~ have several interesting and 
unique properties. Consider for example a pure Yang-Mills theory in D dimensions, which 
has a boson multiplet and a fermion multiplet. Since the longitudinal momentum operator 
is a kinematic operator, the supercharge must be quadratic in the fields, while the 
supercharge Q~, whose square is the Hamiltonian, in general has both quadratic and cubic 
terms in the fields. The dynamics are carried by the cubic terms in the supercharge in the 
sense that a theory with only quadratic terms in Q~ will be a non-interacting theory. The 
supercharges contain an odd number of fermion fields therefore the fermion fields must be 
periodic^. 

We have formulated a number of supersymmetric theories imposing periodic boundary 
conditions on the fields and compared our results with other numerical results. We find 
excellent agreement and we find the SDLCQ converges extremely fast, much faster than 

1 We will not consider the possibility of twisted boundary conditions here, but clearly this is an interesting 
direction to explore. 



1 



any version of standard DLCQ || [TO]. Among the most interesting results we found were 



a number of exactly massless bound states in some theories. These are states that are 
destroyed by one supercharge, Q~, but not the other, Q+, since they have finite momentum. 
These massless states persist at all values of the coupling and it is clear that these states are 
BPS states which saturate the bound \Z\ = M where the central charge Z is zero. It would 
be extremely interesting to find BPS states with non-zero masses numerically. 



2 Formulation of the Theory 

In a light-cone formulation the algebra of the supercharges with a central charge extension 
takes the form 

{Q-,Qp} = PS^ (1) 
= P + S a , 

{Q+,Qp} = Ptffa + ZSaji, 

where P + is the longitudinal momentum, P± is the transverse momentum and P~ is the 
Hamiltonian. Z is the central charge extension and we have suppressed the spinor indices. 
In light-cone quantized field theories one has transverse boost invariance so one can always 
work in a frame where the total transverse momentum is zero, thus P± on all physical states 
can be taken to be zero. 

It is well known that the central charge extension Z can be written as an integral over the 
boundary of the space, and it will therefore vanish if one uses periodic boundary conditions. 
Therefore the last anti-commutator in Eq. (|IJ) always takes the form 

{Qt,Qp} = 0. (2) 

Without a central charge extension the BPS states of a theory will simply be massless states. 
The mass of these states is protected by the BPS symmetry and they will remain massless at 
all couplings. We have seen these states in a recent SDLCQ calculation in 2+1 dimensions 
||. This theory is easily extended to an J\f=2 supersymmetry, but without a central charge 
there is no hope of seeing massive BPS states. A massive BPS state would be a very striking 
feature in the spectrum the theory, since it would have a fixed mass as a function of the 
coupling. 

We propose therefore to extend the definition of the supercharge in order to include anti- 
periodic boundary conditions. This extension breaks the supersymmetry at finite resolution, 
but at infinite resolution it will be restored. We define momentum-shifted analogues of the 
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standard supercharges, Q, i( which carry momentum The Hamiltonian P is defined 

by 

{Q; iQ) Q:i^} = 2v / 2P-^, (3) 

with the normalization of Ref. To test this method in praxi, we study the eigenvalue 
problem 

2P + P~\ip) = M 2 \ip), (4) 

defined by the Hamiltonian, Eq. (|3]), within the context of the supersymmetric theory of 
adjoint fermions in 1 + 1 dimensions. This is one of the simplest possible supersymmet- 
ric theories and all the low energy bound states of this theory are well known |16|. The 
supercharge density is q~(x) 



q (x) = ig2 7 ^Tn 



ip{x)il){x)—il){x) 



(5) 



where ip(x) is a real adjoint fermion field. We expand the field into its modes which will 
carry half-integer momenta, because we impose anti-periodic boundary conditions 
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We define Q.i to be 
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Since we are using a real representation for the fermions it is easy to verify that the Fock 
space matrices Q~i satisfy, 



If all the wavefunctions of the physical states vanish sufficiently fast for small momentum, 
then for large resolution one would expect that the half unit of momentum which the super- 
charges carry is negligible compared to momentum the partons carry. Thus, Q~x — > Q~ in 

2 

the continuum limit K — > oo. However, since we have broken the supersymmetry one might 
face a mass renormalization which is generally required in these two-dimensional theories. 
In spite of this potential difficulty we proceed to calculate the spectrum of this theory, and 
will find that the theory indeed need not be renormalized. 
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3 Numerical Results 



In order to obtain the mass squared eigenvalues M 2 of the theory, we have to solve the 
eigenvalue problem, Eq. (f|). We use the standard large iV c discrete Fock basis. Its states 
are of the form 



\pi,P2, ■■■,Pn) 



Tr 



B\ Pl )B\p 2 )---B\p n )\ |0>, 



where the symmetry factor s counts the number of times the set of momenta (pi, . . . ,p n ) is 
mapped onto itself under cyclic permutations. When using anti-periodic boundary condi- 
tions, the fermion bound states will lie in the sector of half-integer harmonic resolution K 
and boson bound states will have integer K. The action of the operators, Eq. (^), can be 
readily evaluated 



Q + i\Pl,P2, - ■ ■ ,Pn) 



2 7/4 y „ 2J 



Pi±l / I 

r<m) £U+s 



+( 



1 1 
+ — 



vPi-1 Pi 



Pi + Pi 



i-^rj \p,+Pi-i ± i,i, < r"(Pi,p.-i)>} 



where is a permutation of n — k momenta 

Ln\Ph, ■ ■ ■ ,Pi k ) = {Pl+j,P2+j, ■ ■ ■ iPii-l+jiPh+l+j: ■ ■ ■ ) Pi k -l+j) Pik+l+j ) ■ ■ -Pn+i)- 

This compact form allows for an easy computer implementation. Since the supercharges 
change fermions to bosons and vice versa, they change the resolution K by ±|. Thus in this 
discrete representation the supercharge matrix, 

(K\Q- ±2 _\K T \), 

will in general not be a square matrix. The Hamiltonian, however, is constructed as the anti- 
commutator of the of the two momentum-shifted supercharge matrices, Eq. (^). From Eq. (§) 
it follows then that the Hamiltonian is a hermitian matrix. In this matrix representation 
the Hamiltonian at resolution K will receive contributions from intermediate states with 
resolution K ± |. The supercharge matrices in a (fermion, boson) 1 basis have the structure 

f A ± i\ 
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The Hamiltonian for the fermions is therefore, 

Pj erm = A + ,B_ h +A^B +h , (10) 
and the Hamiltonian for the bosons is, 

P^ = B + xA_i+B_iA + i. (11) 

These two matrices are subsequently diagonalized to yield the fermionic the bosonic masses. 
We note that these matrices have in general different dimensions. 

In a first step, we calculated the fermionic and bosonic spectrum up to harmonic resolu- 
tion K=10 by solving the eigenvalue problem, Eq. (^). We then focused on the six lowest 
eigenvalues whose eigenfunctions for resolutions up to K=10 are solely built out of states 
with at most four and five particles for the bosonic and fermionic spectrum, respectively. A 
truncation to four particles allowed us to go to resolution K=25 in the bosonic sector with 
relative ease. In the fermionic sector of the spectrum, the larger dimension of the fermionic 
matrix, Eq. (|TD|), prevented us from going quite as high. 

The resulting bound state masses in units of g 2 N c /n are shown in Fig. [I]. We fitted all 
six curves to the function 



V -r(i\l— + b—. 



We find the continuum masses to be 

r2 



M A aS DLCQ = 22.7, 43.0, 57.7, 58.0, 63.4, 67.0, (12) 
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Figure 2: All the curves correspond to the lowest mass state (in units of g 2 N c /n) of the 
theory. The boxes and triangles represent the lowest anti-symmetric SDLCQ boson and 
fermion eigenstates, respectively. The middle curve stems from a calculation using ordinary 
SDLCQ and the fermions and bosons are exactly degenerate. The lower curve is the standard 
DLCQ result for the lowest mass boson. The upper and lower curves are both fitted to 
M 2 + a'K- 3 / A + b'K- 3 / 2 . 



compared to the values 

M SDLCQ = 23 - 8 > 47 - 9 > 62- 1 ; 62 - 6 > 63 - 8 > 64 - 7 > ( 13 ) 

of ordinary SDLCQ at resolution K — 8. While the agreement is far from perfect, it is clear 
that for this model the anti-periodic SDLCQ method works, in the sense that it correctly 
reproduces the spectrum of the theory. It is fair to say, however, that the exact continuum 
values are quite sensitive to the fitting function. 

In Ref. || we compared the lowest eigenvalue as calculated in SDLCQ and standard 

DLCQ as a function of the resolution K. We found that the convergence of standard DLCQ 

is slow and clearly nonlinear as a function of 1/K. Unfortunately, we see from Fig. |2] 

that the convergence for anti-periodic SDLCQ is, roughly speaking, about as bad as the 

standard Hamiltonian DLCQ methodF]. In Fig. ^| we show the lowest fermion bound state 

mass and the lowest boson bound state mass as a function of 1/K for the anti-periodic 

SDLCQ approximation. While the fermions and bosons appear at different resolutions, 

they can be fitted by the same curve for large enough K. Thus the supersymmetry of the 

2 For comparison, we used the same fitting function that was used in the DLCQ calculations of Ref. || 
for our data. 
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spectrum appears to be recovered by the approximation already at low harmonic resolution. 
Unfortunately, although the potential perils of supersymmetry breaking seem to be absent 
already at K — 10 for the lowest eigenstates, the convergence remains unpleasantly slow. In 
principle, however, the method appears to be working, and we anticipate room for numerical 
improvements. 



4 Conclusions 



In the present note we introduced a novel version of SDLCQ that allows for the use of 
anti-periodic boundary conditions for all fields. Consequently, the formalism can support 
for the first time boundary integral contributions. We provided evidence that the approx- 
imation converges to the standard supersymmetric results. The new approach inherits the 
very advantageous feature of absence of mass renormalization from SDLCQ. Unfortunately 
however, the approximation does not enjoy the same rapid convergence as SDLCQ. The 
lack of rapid convergence is a severe problem, since supersymmetric theories with non-trivial 
central charges (N > 2) have many species of particles and therefore a very large Fock space. 
Thus going to sufficiently high resolution to obtain accurate results for the spectrum will 
be quite difficult. Of course, the standard DLCQ approach can be used with anti-periodic 
boundary conditions. But, as we saw in Fig. |2|, its convergence is no better than that of the 
anti-periodic SDLCQ approximation. Moreover, the standard DLCQ approach will certainly 
have renormalization problems in higher dimensions. 
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